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SUMMARY 



The present report deals with the variation of a tur- 
bulent velocity profile in flow from rough to smooth trail 
and vice versa. Expressions obtained for the shear-stress 
distribution with respect to the distance from the point 
of junction of the different roughnesses and from the wall 
distance, are utilized to ascertain the developing veloc- 
ity distributions. Under simplified assumptions, the use 
of these formulas renders possible the integration of the 
motion equations for the shear stress. This calculation 
is carried out and compared with the experiments. Despite 
the fact that the assumptions in this particular case do 
not prove to be wholly correct, comparatively good agree- 
ment is achieved in the most important region. 



INTRODUCTION 



The purpose of this report is to study the variation 
of a turbulent velocity profile established in a channel 
of known surface roughness on transition to a different 
surface roughness by experiments, and to develop the for- 
mulas necessary for the calculation. The phenomena accom- 
panying fully established flow along a boundary of known 
roughness may be considered as being' understood now (ref- 
erence l). Expressions have been obtained for smooth and 



"Umform.ung eines turbulcnten Geschwindigkeitsprof iles. " 
Zeitschrift fur angowandte Mathematik und Mcchanik, 
vol. 19, no, ?, April 1939. pp, 87.-.100, 
This article forms the second part of the thesis, un- 
der the direction of Professor Prandtl, entitled.: 
"Studien zum Rauhigkei t sproblem . " The first part, en- 
titled : "St romung .hint or einom. o'inBclnon - Rauhigkei t s- 
element," a^opeared in' Ing . -Arch'-*, vol. IX, no. '5, 1938, 
p. 343, 
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rou^;h tubes and channels for the velocity di st rihiit i on , 
resistance, exchange in impulse, and tho mixing path. 
If the surface rou^hnoss undergoes a- change , a new veloc- 
ity profile has to he formed, and this transition region 
has, so far, heen investigated only theoretically, under 
certain simplifying assumptions (reference 2)» The au- 
thor ha.s carried out experiments over this region and 
utilized the results to esta'blish a new method for com- 
puting the velocity distribution over the transition sec- 
tion of a channel when changing from smooth to rough, and 
vice versa, provided the shear stress at the wall under 
conditions of fully established flow, is known. The pre- 
dictions have "been confirmed "by further experiments. 

EXPERIMENTAL. SET-UP 



Experiment s .were carried out on a "blower whose nozzle 
emptied into a rectangular channel 60 cm wide and 20 cm" 
high. The channel itself consisted of airtight wooden * 
"boxes, holted toge'ther with felt strips in between, insur- 
ing a completely closed channel length which could be 
lengthened or shortened as necessity arose, Tho end of 
the channel formed tho experiment chamber. Glass windows 
on the side assured the most accurate setting of the pitot 
tubes. In virtue of the reflection on the smooth wall, an 
accuracy up to l/lO mm was readily obtainable for a known 
point of contact of the survey tube with the wall and 
hence, of the \vall distance. The experiment chamber itself 
was so arranged that the survey tube could be moved across 
the entire length of the box but not crosswise to the direc- 
tion of flow (fig, l). The box was .lt5 m long. 

Above it was the guide rail of an optical bench, to 
which the streamlined tube. was fastened. The survey tubes 
were fitted with threads and laterally screwed to the 
streamlined tube. The recorded pressure was transmitted 
by means of a valve rubber to a little tube inside the 
streamlined tube and was removable from above. The slot 
in the upper wall'was covered during the test by well- 
fitting wooden strips. 

Since a two-dimensional problem w.as involved, the 
static pressure could be recorded with. a disk type of sur- 
vey tube originally employed by Motzfeld; (reference 3), 
consisting of a lens-shaped disk p-f 8 mm 'diame t er , 1. mm 
maximum thickness, with 0*6 mm diameter orifices in the 
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cent er. of . the" "curved stirfaces. . The . survey . tu^be should "be 
so introduced that the mean* flow direction is always par- 
allel to the.tuhe. The ' di recti onal ' unsusceptibili ty was 
approximately ±3^ and v/as at any time ohtainatle. A call- 
"brat i on . was made in the free stream, hack of the fan , in 
the "sound" flow, with the static pressure set at zero. 
The coefficients were almost i dent i cal wi th the employed 
survey tubes, averaging around 0 =1.14. This value was 
lar-?ely independent of the Sey^o^^s numher. With the no- 
tation : 

Pg total pressure 
^st st^-tic pressure 

P' '"static pressure recorded 
^ hy the survey tuhe 

q[ dynamic pressure 

We have : 

Pp- " Pst' 1 

Pst = - 1 = Ps - ---p = p - 1) Pst'^l 

which, with g = 1.14, 5;ives the sts^tic pressure at 

^st ^ 1.14 

All pressures were recorded with Prandtl manometer 
relative to the pre'-ssure in the experiment chamher. Strips 
0,7 cm hiii^h and 1.0 cm wide, nailed 15 cm apart on the 
bottom side of the channel, simulated the wall roughness. 
The rough length amounted to 5m. 

PRELIMINARY TESTS 



Before proceeding to the actual measurements, it was 
necessary to ascertain v/hother a fully established pro- 
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file existed at the end of the rough^ length. Plotting the 
velocity; a^jainst the logarithm of the wall distance, ^*ave 
a strai'-?ht line which, acco^din^5 to the logarithmic law of 
velocity distribution, is an indication of fully developed 
flow, as was to he expected on the "basis of Nikuradse's 
earlier pipe flow experiments, which established the ter- 
mination of mixing after about 40 pipe radii, Kirsten 
(reference 4) claims the presence of the final velocity 
distribution at 20 pipe radii already. With a hydraulic 
radius of 



rv = 



2 channel section 
with circumference 



= 15 cm 



in our case we obtain a distance of 3!5,3 radii available 
as entrance which, considering the marked roughness, should 
be sufficient for full development of the profile. 

The measurement of the pressure drop was accxjmpanied 
by a slight pressure jump behind the roughnesses, on 
transition to the smooth wall, caused by the abrupt widen- 
ing of the channel. A smooth, carefully alined plywood 
board, high enough to assure progressive transition of the 
pressure, was nailed at the smooth part of the length. If 
the board was too thick, naturally the opposite occurred: 
a pressure jump to too small values. After various trials, 
1.0 cm was chosen as a practical height. 

The static pressure drop records were made in the 
middle of the channel. The effect of v, that is, of the 
velocity perpendicular to the wall, on the pressure read- 
ing., was below the accuracy limit; hence, could be disre- 
garded. Tor, with P^t^ as the additional pressure indi- 
cated by V, our calibration above affords: 



st. 



- V 



st 



» = - 0 



p 

,14 I V 



^st 



st 



0.14 



u 



a 



Assuming roughly that v/u in the center acts as* 1:100 - 
this value is not reached in our case - we have: 



st. 



- P 



st 



Pst ' - Pst 



1 

10000 



The omission is therefore fully justified. 
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The slight i r ref^ulari ty of the velocity distribution 
across the width of the channel, was of no import in these 
tests as the iaeasureinent,s were made in a vertical plane • 



TEST PROCEDUHE 



Having to do with a two-dimensional problem, the x-y 
coordinates are formed hy the perpondi cular plane of sym- 
metry of the channel in the flow direction, its zero point 
being situated in the junction point of the two roughness- 
es. The por.itivG x axis points in the direction of the 
principal flov;; the positive y axis, perpendicularly up- 
ward. 

The flow was studied when changing from smooth to 
rough wall and vice versa. However, for reasons evolving 
from the experiments, only the variation in the profile 
from rough to smooth, is to "be analyzed in detail. Veloc- 
ity measurements were made at x = 2, 7, 15, 20, 40, 70, 
100, 150, 215, and 290 cm (fig. 2), although only part of 
the obtained curves are shown, for the sake of clarity. 
The variation in the profile is such that in the lower re- 
gions with increasing x, an increase - and in the middle 
a decrease - in velocity occurs, whereby the smooth pro- 
file evolves gradually. This transition is numerically 
treated. 



ATTEMPTED SOLUTION 



With T = apparent shear stress, the equation of mo- 
tion reads: 

u 5il + V Sii - ^ 1 ^ 1 3 T 
ox 3y " p bx P By 



. With the continuity equation ~ + = 0. we obtain, 
conformable to a variation by Prandtl (reference 5): 

. U -l^ 4- V ■= - i ^2 + i = f (y) 
dy p 8x p- dy ' 

where, for a fixed x, the right-hand side is put equal 
to a function of y; or 
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We integrate from 0 to y , since v = 0 for y = 0, 



u 



y 

--^^■5 dy ■ and hence, v = - u / —^^g dy 
u(y) ./ u(y) 



whence differentiation :?ives the .velocity increase: 



- ^ = ^ Bx ^ ay C^y 

0 . 

Starting therefore from a certain velocity profile., if -g^ 

(which i^n first approximation can be put constant because 
of the small curvature of the streamlines over the channel 

secti'on) and - are known, ^ can "be computed, wheVe- 
from follows a new u for (z + Ax); 

The method applied recurrently ^ives, therefore, the 
velocity distributions at the various distances- For the 
shear stress, the formula 



T a I 3u 

P " ^ I ay 



^ (reference 6) 



was applied. It could be computed under rou^h assumptions 
of mixin? path I. Reliable mixing-path distributions in a 
channel are unknown. Hence, the use of the simple formula 

I = 0.4 y which. for y = ii five's,,. I = 0.1 h (fi?. 

h 2 

3); h = channel height. 

The calculation was then attempted with -.thi s • assump- 
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tion* The velocity u. of the initial - prof ilo supplied 
the measurement, had to he defined ^graphically Vy 

differentiation.' The result V7as the shear stress and 



9t 

after another differentiation -5-. The inte-?ral 




0 



dy was graphically evaluated. Reiterated differentiation 



of u / dy finally gave the desired velocity dif- 

' ^ 0 

f erence for the selected path element A 

Notwithstanding the various applications of the usu- 
ally not too-acciirat" e differentiation and /graphical inte- 
gration, we still "believed that with the care use*d, the 
results would agree with the actual conditions; hut it* was 
otherwise. A coraparison ^vith the experimentally ascer- 
tained velocities manifested, to "be sure, close agreement 
for greater distances from the wall, "but also great de- 
partures in the neigh'borhood of the wall. The advance of 
the rough profile toward the smooth, up to y values of 
several centimeters, took place much faster than the cal- 
culation stipulated* The reason could only be ascribed to 
the erroneous mixing-path assumptions. 

In the attempted explanation of these conditions, the 
opposite process was essayed - that is, the shear stress 
was computed from the recorded velocity profiles and their 
differences, and the extent to which the aoove functional 
relation of mixing path and y was true or false, and 
was independent of x, was checked "by means of the formu- 
la - = l""! -^i ~ , The equations of motion and continu- 

ity served as basis. 

The result definitely ostabii.shed the , reason for the 
discrepancies between calculation and experiment, a^nd made 
the continuation of the above method appear little prom- 
ising. The rise of the mixing path with' the distance from 
the wall was much more pronounced than we had stipulated 
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in our forraula, according to Nikuradse^s experiments on 
pipe friction. But. at greater x values (150-200 cm)..- 
that is, on approaching the fully esta'blished smooth pro- 
file of Nikuradse 's-.distri'bution^ the mixing-path;, distri- 
hution was approximately similar (it was a little "below). 

Incidentally, it may he mentioned that, according to 
'Prandtl, the mixing path on approaching the wall ~t ends to- 
ward a certain limiting value, which corresponds to the 
degree of roughness. This value, "being small, was dis- 
regarded in our calculations -To illustrate: Consider our 
utilized rough wall; writing, in conformity with the uni- 
versal law of velocity di st rl"bution for rough walls (Nik- 
uradse, reference l), 

8.5 + 5.75 log 7^ = — = 0 

gives the wall distance y^ , with velocity = zero, accord- 
ing to our law. v* is the shear stress and k^ the 

r s 

equivalent grain size, which will he determined later. .As- 
suming its value for the present, gives: 

1„ Li , log ^- = - y„ = 0.175 cm, V„. = 0.068 cm 

s 

This value 1^ would not have heen a'ble to change the mix- 
ing path distribution very much. On the smooth channel 
wall used it was, in fact, several percent lower, so that 
it could not he considered at all* The wide discrepancies 
therefore remained unexplained; hence the solution of the 
posed problem hy this method must he called unsuccessful. 



FORMULATION OF AN EMPIRICAL SHEAR STRESS 
DISTRIBUTION FORMULA 

a) Determination of wall shear 'stresses and closer ' 
characterization of wall roughness^ - Before g:oing on with 
the new- method, the utilized roughness and it$ shear stress- 
es are scrutinized somewhat more' clo seiy. * , 

T is the shear stress on the smooth wall ; that 
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on the rough wall;. T :. was o^btained once "by pressure drop, 

at ~ = 4,90 X ib^.cja^/s^. But since this value seemed 
too uncertain, the determination of "by the universal 

law of velocity di st ri'bution : 

u g 
= 5.5 + 5 . 75 log 

. -J 

served as check. Logarithmic plotting of the wall distance 
y / against the velocity u affords a straight line* 

u = ffig + ng log y 



from whose slope r^^ follows direct 

ng = 5.75 v*^. v*^ = = 



It gave: 



T 

= 5.05 X 10^ cm^s^ 



The factor 5.75 was considered safe. The two ohtained val- 
ues agreed fairly closely. The additive constant of 4.32 
in our universal' velocity distribution law differing from 
Nikuradse's 5.5, is at t rihutahle to the fact that the wall 
is not perfectly smooth. Thus the universal law for our 
case • reads : ■ ■ 

V* y 

— rr =: 4.32 + 5.75 log — ~ 

^*g ^ 

The slope of the . straight s is the same as Nikuradse's. 

With a = side length, h = height .of the channel 
section, T^^/ p is now computed from the force equilihrium 

"between wall shear stress and pressure gradi exit : 
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Analogous to the determination of. Tg/p, Tj./p can now 
he ohtained for checking the last-ohtained value. For a 
rou^h wall, the universal velocity distribution law takes 
the form - 

— = A + 5,75 log* -21./ u = m- + n^ lo^ y 
r s 

which, plotted as previously, loaves 

= v^^ =: ^ = 155 cm/s. ~ = 24.1 x 10 cm /s 
Even in this case, the values are fairly a^reeahle. 

Tp »j 2/2 

With the approximate average value -r = 23.0 X 10 cm / s 

r 

were computed* 

It should he noted that the roughness of the hoard 
used for the cross-sectional contraction, was slightly dif- 
ferent from the channel wall* This explains the differ- 
ence between the above-defined wall shear stress, applica- 
ble for the usual channel wall, and that later obtained as 



1 n 9Tdy 

constant after integration of p / """^y""* *which represents 

J 0 

T 

the shear stress of the board.' This value ranged from - = 

3 X 10^ to 4 X 10^ cmVs^. Subsequently, 3.4 x 10^ was 

used in the calculation as it involved flow conditions 
over the "smooth" board. 

At this time, a remark about the shear stress on the 
smooth wall when facing the rough wall , may not be araiss. 
From the smallness of the discrepancies appearing in the 
measurements, it may be concluded that in this case the 
wall shear stress is -greater than when the whole channel 
is smooth. The rour?hness effect extends across the entire 
channel . width up to the opposit e wall. 

For comparing our roughness with that of others, we 
effected a reduction to Nikuradse's grain size, although 
this should not be looked upon • as generally the best ref- 



NACA Technical Momorandum No. 951 



11 



erence roughness. One of Schli chting • s ?rain sizes (ref- 
erence 7) mi^ht perhaps be more appropriate since in . 
Nikuradse's rou.?hness hecause of the adhesion of* the 
sand grains V7ith lacQuor - the conditions "between the 
grains cannot he accurately simulated. 

The equivalent ^rain size ^ives the grain si%e of 
Hikuradse's sand roughness, having the same resistance as 
the roughness used; it is -indicated .hy kg. The universal 

velocity distribution law jTor rough wall for fully estao- 
lished flow, reads: 

u y 

= A -f 5.75 log f 

V* k 
r 

A is characteristic of any roughness; it is a roughness 
function.. For sand roughness, the value is Ag = 8,48 . 
Accordingly, 

-™'= 8.48 + 5.75 log 7^ 
V* k^ * 



and, comparison with the preceding relation: 

8.48 + 5.75 log = A + 5.75 log f 

kg ic 

or, after combining: 

ks 

5.75 log — = 8.48 - A 
k 

Sorailogarithmi c plotting gives, in our case, A = 3.45. 
Hence , 

k k 
log ™ = 0.876. = 7.52 . 

k . k 

which, for the k 0.7 cm height of roughness used, 
amounts to • 

kg = 5.26 cm . ' 

The size of the equivalent sand roughness is therefore 
7^ times greater than that of the used roughne ss . 
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Further cal culat i on . demands an exact definition of 
the wall distance. Schlichting defined it as equal to the 
distance of a fictitious , wall ' suhstituting for the- rou^'h 
wall, havin<? the same fluid volume. Thi's.v/as found:- to -he 
impractical. The starting point formed the lo^jarithmic 
veloci tj^-di st rihut ion law : 

y 

u = 5.75 log ™ 

•^0 

y 

Ic?; — was then plotted a^-^ainst u^ while .y^ was var- 

^0 

ied until the extension of the obtained straight lines 
passed throiigh zero. The points adjacent to the wall were 
omitted, since the law does not apply there. The thus- 
obtained point y^ represents the heif^ht at which the ve- 
locity is zero, and consequently, forms the i^ero point 'of 
the new coordinate system. y^ was equal to 0.17 cm. 
This agrees with Nikuradse's definition, according to 

which y^ = In our case k = 0.175 k. 

The velocity measurements were made at an average air 
speed of around 14 m/s. The velocity (16.40 m/s) recorded 
in the middle of the channel served to effect the nondimen 
sionality, 

b) Mathematical shear stress di st rihution. - Reverting 
to our original prohlem - that is, the calculation of the 
variation in velocity profile - it was necessary to oh- 
tain some kindsof expressions. A ray of light in this di- 
rection v/as indicated during the plotting of the shear 
stresses computed from the experiment. The shear stress 
"by fully estahlished flow is, as known, linearly distrihut 
ed across the section. It is zero at maximum velocity and 
increases with approach to the wall up to the correspond- 
ing wall shear stress. The shear stresses prevailing at 
different distances from the roughness transition point, 
were all "between the linearly rough and smooth distri^hu- 
tion, working from the rou?h with increawsing x oyer to 
the smooth (fig-. 4). The new wall shear stress was imme- 
diately availahle. Its effect spread consistently upward. 
Hence, it was presumed that the differences "between the 
original rectilinear shear stress di st rihiition and that 
existing at certain distances, might oe 'approximately ex- 
pressible hy an exponent ial. relationship . he.tween^ -y and 
Xt Comparisons of various formulas with the obtained 
shear stress distribution, manifested the following as the 
most favorable : 
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T(x.y) =CTj. - ( Tj. - T^) f(Tl)]^~-2f- 



(1) 



f(Tl) = a e""^^ with T) = \ 



(2) 



y = h indicates the ordinate of the maximum velocity. 
For fully established smooth and rou^h flow, h assumes, 
of course, different values - a fact which was not taken 
into consideration, as its role is subordinate in our anal- 
ysis, A certain compensation mi^sht he achieved hy resort- 
ing to the arithmetic mean of the hei£;ht as oasis of the 
calculation. Our choice on transition from rough to smooth 
as heli^ht h is the distance of the velocity maximum on 
the rou?h profile from the wall. 

During the closer investigation of our previous sur- 
mise f(T1) was first plotted from equation (l) against y 
for different x (fig"l 5): 



f(T1) = 1 + 



T(x,y) h 



- T h . y 



A certain value f(T|) = a e"*"*^^/^^ = constant defines a 



certain vdlue Tl =• = constant 

. x^; 



Therefore , 



" In y = In const + m In x 

After logarithmic plotting of y and x, followed "by 
the drawing, of straight lines through the points obtained 
for f(Tl) = const - the errors encountered "being compensated 
for 8.S much as possible - their slope ra indicates the ex- 
ponent of- -X (fig.. 6). The gratifying feature of our for- 
mula (above) is that, with it the average slope of the 
straight lines was quite constant; i.e., approximately 
equal to 1, in consequence of which. f{T\) assumed a form 
easily amenable to calculation. T\ forms the intersection 
on the ordinate obtained ,by extension of the straight line 
placed as closely as possible through the points with the 
averaged slope 1 as far as the intereect ion with the ordi- 
nate axi s . 



There "remained the determination'' of ^- a and b: 
f(Tl) = a e ^, log„[f(.Tl)l = log a - b T] 
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SemiloJ^arithmic plotting then " gave ^ *a' arid h , ' For a we 
obtained about 1, as it should "be * if T^^^^ " g 

small y; for h, we obtained 11-. 16 (fig. 7). 

Then, the general relationship of x and y for T 
reads : 

U.y) - (Tr-T,) e-"-"s]^-5-^ . (3) 



This is an empirically established formula for the change 
in shear stress distribution on transition from rou^h to 
smooth wall, 

Kow it remained to be proved, to \vhat extent a -gener- 
alization of this formula was permissible to other cases.. 
For this purpose, the in vesti #?;at ions . were made by exactly 
opposite conditions, as in flow from 'smooth to rou^h wall. 
The length intended as entrance obtained a bottom support 
of 1 cm thickness, of the same boards used in the first 
test and, consequently, with exactly the same roughness, ■ 
The entrance length of 5 m was considered sufficient for 
establishing the profile. 

Next comes a rough length of about 3 m. Great care 
was used in fastening the board on the bottom, so as to 
obtain a, perfect plane and keep the cross-sectional 
changes within negligible limits. The smooth piece was 
followed by the strips spaced at the given distances. The 
measurements wore made in the exact center between strips, 
x= 7,5, 59r5; 71.5, 103.5, 140. 204, and 290 cm (fig. 8). 
Jot the rest, the calculation was as before. Figure 9 
shows the shear stress distributions. Figure 10 contains 
the representation of fCH). The constant factor in ex- 

-b^ 

ponents of f(Tl) = a e ^ disclosed a slight difference 
from the previous value (fig, ll).. But this error was too 
insignificant to entice us into attempting to express both 
processes by the same formula (fig, 7). The mean value of 
b represents the above value 11.16. a again followed at 
around a = 1, and this value 1 was used thereafter. 

Our. shear stress distribution formula for flow to 
rough wall, has thus the same form as for flow from rough 
to smooth wall, except for the exchange of and T : 

r -^11., 16 - . 
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The same type of . shear stress yariat ion*. testifies to the 
same type of turbulence. ' 

The extent to which this formula actuall;/ . repre sent s 
the experimentally achieved distribution, is shown in fig- 
ures 4 and 9. In the usually most important region of 
small X, whore transformation advances most rapidly, the 
assimilation is almost complete ; subsequent discrepancies 
are of much less influence. Concerning the small differ- 
ences of the experinentally defined shear stress distribu- 
tion in wall proximity, it nay be stated that they ar$ prin- 
cipally due to the. fact that the integration constant s ' re- 



suiting from the integration of ^' p 8y ^representing 

'Jo 

the shear stresses on the wall, did not produce exactly 
identical values. The^^e differences are caused by inaccu- 
racy of moasiirement and evaluation. In this re^^ion the 
formula seems to reproduce the conditions much better, 
since the sa,me wall shea.r- stress must prevail. 

In retrospect, it m.ay be stated that the above formu- 
la reproduces the conditions fairly well and can be used 
as a basis for oompxitin^ similar problems. In general, the 
method of differentiation will have to bo rosorted to - 
i.e., start from a known, fully established profile and 
define du/3x, ..exactly as in the first attempted solu- 
tion - this time with the. new shear stress formula - ob- 

tainin^, according to Uo = u, + I :r-- "^i Ax a new -orofile. 
: ^ Vox /i ' 

Graphical treatment produces quickest results. The x in- 
tervals need not be chosen so very small in order to 
achieve a ?;Dod reproduction of the actual curve. In our 
case, sections of from 10 to 20 cm directly behind the 
junction point should be chosen, which could even be in- 
creased as the spacing is increased. 



SOLUTION OF MOTION EqUATION WITH SIMPLIFIED ASSUMPTIONS 

In the following, it is attempted to integrate the 
m-otion equation .under simplified assumptions by introduc- 
tion of the shear stress formula. The equation of ..mot ion 
reads: 
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ax 3y p.6x. P ay . ' 

Our previovisly obtained formula serves as "basis: 

' ■ • y 

T = ^-~r[T, .+ (T, - T,)e" ^3 

a = 11.16. and • denote the equililDrium wall shear 

stress for x < 0 and x + cd. 

Now the first simplifying assumption stipulates: 
I Tj^ " ^i-l ^"^1 • Hence, vre can write 

u(x,y) = u*(y) + u'(x,y), v(x,y) v«(x,y) 

with u' and <. U. u^(y) is the velocity profile 
over the roughness 1. With the stream function 

Bu ' da * 

while disregarding the small terms v.' and u' g~ , 

we ohtain approximately: 

^ .Alt . ^it ^ ^ll2 ^ l-ll (7) 

* Bydx ay Bx . P 9x p Sy 

With equation (6), the right-hand side assumes the form: 

a y .a^ 
p 5x phL^i + ^T2-T3.; e J ^^Q-^i) ^ ^ 



Putting 



p Sx P By P 3x ph 



that is, f(x) = 0 for x< 0, the right-hand side of 
equation (7) gives: 

* • 9". y ' ' ' ; . • 
- (i * a (T,.T,) o' ^ - - - - ■■■■ 
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For mbdiiua x and not unusually groat; • a(h - y) » x 

for X = ay, tho expression is 

:- x .+ a(h' - y) ' r: a ■ ■ ' ■ ■ 
•"h X X 

which is to sorve as mean value in tho subsequent calcu- 
lation. This comprises the. most usually interesting 
region (x = 11 y ) • "The velocity* distribution u^ it- 
self is to be represented by an intervening straight lineJ 
u# = A + By. The selection of this straight line should 
be such as to assure .its close approximation to tho veloc- 
ity profile in the pertinent region. 

Substituting v for - and posing for abbrevi- 

c)x 

T — T 

at ion for tho constant value, ~ ^ = C, equation (7) 



gives 



- (A + By)|^ + 3v = f (x) - J C e-^^'x (g) 



This is a linear inhonogenoous differential equation of the 
first order, v/hich is solved by the method of variation of 
constants. The solution leaves: 



y f (x) C e ' 



- V = (A + By) P ^ dy (9) 

J (A + By)^ 

0 

The integration constant is zo.ro , since on the v/all, that 
is, for y = 0, v oust equal zero. The condition of 
continuity on the other wall is net by means of f(x)» 

For the further treatnent of the above integral 
(A +''"By)^ '' developed by binonial expansion, whereby only 

the first two terns are considered. Thus. 



£ \ . A / 



(A + By.) 2 A^ ^ .A 

and therefore 

y 

A^ 



[f(x) - f C y/.^ j (l - ^) dy (10) 



18 



UACA Tochnical Momorandun Ho, 951 



ThG integration affords 



A2 



f(x)(y-|y=).Oo-?(l-fy-f t>3f t-oj 



(11) 



whence the continuity equations give the velocity u at 



u = 



- r ^ dx + u. 



In conformity v/ith the initial condition: x =r 0 , u =: U, 
the integration constant inust he eq\ial U. Consequently, 
-V* must "be differentiated again v/ith respect to y and 
then integrated v/ith respect to x, which finally leads to 
the integral logarithm and, after intermediary calculations: 



u 



B C 



LA a 



CO 

- r 



5 3' 
A3 



1 P 

•y^ i / f (x)dx + X e ^ 



[ 



B C C 



(!-)] 



■s dz I — 



Ca 3 B^'C 



C 2I 
- a y + U:, 



(12) 



2 = v/as chosen as a substitute. Since the integral 

logarithm produced is tabulated, the prohlem may be consid- 
ered solved* 0?he integral logarithm is usually expressed 
(see Jahnke-Emde) by: 



c» 



/ 



r i_i dt =-E.(-x) 



(13) 



X 



v/hich, in o'ur. case, gives: 
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00 




z ■ 



dz ■■■ 



(14) 



In this manner it is possible to compute, "by known 
wall shear stresses, every profile formed on transition 
to Tg , ' starting from a velocity profile formed over T^, 
ifc proceeded from the premise that | T3 - | <: • 
It should, however, he eciually of interest to apply this 
calculation once for our experimentally explored case, 
notv/iths tanding that the assumptions are not completely 
satisfied and discrepancies therefore may he expected^ 

The initial profile vras that shov/n in figure 12. 
The intervening straight line serving as a "basis of the 
calculation is chosen so as to assure comparatively good 
agreement in the lower region of the velocity profile. 
Surprisingly, there is a fairly clpse. accord with the 
measured velocity profiles for small y. At greater v/all 
distances the differences arc cons iderable. Above all, 
the overlapping of two consecutive profiles docs not ap- 
pear* The effected omissions are evidently too great. 
The principal reason probably lies v/ith the inadmis s ibil- 



Even the approach to 1 causes serious falsification of the 
rosult through exclusive consideration of the first tv/o 
terms. Aside from that, there is also the ommission in the 
friction terms. 

More accurate results, to be sure at greater expense 
of time and labor, are obtained if the complete friction 
fdrnula is allowed for in the integration of the above 



fore, introduced in the calculation through the .first two 
terms in the binomial expansion. This omission- is inevi- 
table if the calculation is to be at' all feasible. The 
solution is : 




for great dis- 



¥ < 1. 



differential equation and 



1 



- alone is, as be- 




NACA Technical Menorandun Uo,. 951 



u =• 



3_e 

A2 



U a J LA 




+ X e 



-a y/x r B g _ 

L A^ ■ 



00 « 



ayz 



dz 



' C g 
L A 



CO 



3 C ^ 3 



+ U 4- 



g-ayz ^ 2 C ay ECay® 

/ da *- — 

L h A h. A^ 



3h A' 



- -^f [ 3^ C /a 2 ^ 7 x^y ^ 2 x^' \ C 1 2B^C_ ^3 
- Lh A= ^3 " 3 a Sa-y h.A J 3 h A^' 

(15) 

Pollov;ing this, the velocity for x = 40, y = 1 v/a,s 
computed as a check on the c£Llculation. The difference 
from the result dDtained v/ith equation (l2) was very small 
as the effect of the additive terms of the last equation 
does not "become noticea*ble before x ^ ay. 5^or the region 
x = ay, equation (l2) v;e.s quite accur?vt elj^* applica-hlc. 
Grreator y vs.lues are, of course, again excluded, "because 
the binomial expansion is then no longer admissible. 

In order, therefore, to include this region of greater 
y values as well, it vj-as attempted to replace the earlier 
oblique intervening straight line of the velocity distribu- 
tion by a distribution U = const, that is, by a vertical 
straight line,ana the choice fell to the average velocity 
u = 1390 cm/s. Ihe degree of the differential equation was 
lov/ered hereby. 

From equation (7) follov/s: 

1 3t) 1 9t . . . /ay 
dydx p Ox p 9y ^xh 



X h/ 

(16) 



V/g find 
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Therefore 



6x 



= 1 r [f(x) .(^ -^-i)c o--^/^]dy 



With the aid of the continuity, v/c ottain for u: 



(17) 

whore donotos the initial profile and u' the addi- 

tional velocity. 

vrith — = z chosen as suhst itution, the integration 
X 

e-^^"" C X -ex y/x 



dz - r—r o 



z Ah 
i/k (18) 

The vclccitj- profiles v/ere computed for x = 20, 40, 
70, 150, and 290 cm distance, A comparison of these re- 
sults v/ith the experimental velocity distributions surpris- 
ingly disclosed a relatively good agreement, especially in 
the middle of the channel (fig. 13). The- overlapping of 
the theoretical profiles already occurs at smaller y val- 
ues, so that here the departure is somcv;hat greater. In 
direct proximity of the v/all the calculation is inapplicable 
becp-use the integral logarithm 

-K,.(,^) = -1,(0) 

becomes infinite. 



To achieve a closer approach to reality, the profile 
could nov; be replaced in the lov/er region by an oblique, 
in the middle by a vertical straight line. A stepped dis- 
tribution v/ould also be fea.siblc, special consicrations 
being then necessary on the points of discontinuity of the 
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velocities* A comx/arison of such results with tho oxpori- 
ment vrould surely be" intexos ting. Y.o.t another profitable 
study should bo the extent to v;hich tho shear stress for- 
mula, obta.ined a,bovc holds true in a comparison with the 
actual processes under entirely dissimilar tyjoes of rough 
ncss* 

Translation by J. Vanier, 

National Advisory Committee 
for Aeronautics* 
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Figure 1.- View of teet cbaxmel. 
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Figure 4.*> Shear stress distribution on 

transition from rough to smooth wall* 



Figure 5.- 
f('7) plotted 
against wall 
distance y 
for different 
parameters x. 
(from rough 
to smooth) 
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Figure Velocity profiles on transition from 

rough to smooth wall. x=distance 
from last roughness. 
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Figure 3.-- Sketch 
of assumed mixing 
path distrihutlon. 



Figure 6.- Log y plotted against log x for different 
f(^) a-e"by/x"const. for defining exponent 
m from the slope of the curves, (from rough to smooth) 
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Fisrure 10.- f (ti) plotted against 

wall distance y for 
different parameters x. (smooth 
to roufi:h) 
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Figure 11.- Determination of 

exponent m for 
transition from smooth to 
rough. 
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figure 8«- Velocity dietrilmtions on 

treuisitionB from smooth to rotigh« 
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Figwe 9.- Shear etresa distrihutibn on transitions 
from smooth to roiigh. 
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Figure 12.- Comparison of calculated and measured 

velocity distributions on the hasis 
of the oblique intervening straight line* 
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Figure 13««- Comparison of computed and measured 

velocity distributions » the calculations 
being based on a velocity distribution 5* =0.85 
constant over the section* 
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